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based on acoustic ray tracing theory as applied to linearly 
segmented sound velocity profiles (SVP) . The models were 
developed into three calculator programs, two for CZ range 
predictions under different source and receiver depth con- 
ditions and one for CZ gain and transmission loss (TL) pre- 
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I. THE NEED FOR GOOD CZ PREDICTIONS 



A. THE PROBLEM OF THE AIRBORNE ASW UNIT 

In most of the Pacific Ocean and much of the Atlantic 
Ocean and Mediterranean Sea, convergence zone (CZ) conditions 
exist a majority of the time, and they provide passive acous- 
tic sensors with an important means of detecting sounds emit- 
ted from submarine targets . In some areas the CZ regions 
are the most important contact regions considered in planning 
acoustic searches. Obtaining accurate predictions of CZ 
sound propagation is therefore vital to the success of acous- 
tic sensor tactical planners. 

Currently, the primary source of acoustic predictions for 
U. S. Navy units is the Fleet Numerical Weather Center, Mon- 
terey, California. Propagation loss profiles for four stan- 
dard frequencies and three source and receiver depth 
combinations are normally provided in the ASW Range Predic- 
tion System (ASRAPS) to air ASW units when requested. The 
profiles, showing transmission loss (TL) versus distance, 
are generated on a large digital computer which uses the 
Fast Asymptotic Coherent Transmission (FACT) model. This 
model uses as inputs a linearly segmented sound velocity 
profile (SVP) , source and receiver depths, and frequencies 
of interest. The SVP used can be specified by the user or 
can come from information stored at FNWC in the form of 
historical data. This stored data is updated by bathythermo- 
graph CBT) reports through a complex weighting scheme as the 
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reports are received. Without going into further detail, it 
can be stated that predictions produced are only as good as 
the data and the computer model used, and only as timely as 
communications allow. 

When entering a search area, a problem often arises con- 
cerning the TL profiles obtained from FNWC. Upon taking a 
BT measurement, the unit often finds the BT profile used to 
generate the acoustic predictions does not agree with the 
actual BT conditions in the area. If this situation occurs, 
the unit tends to lose confidence in the accuracy of the 
predictions and tactical effectiveness is felt to be dimin- 
ished by lack of good information. The objective of this 
study was, therefore, to investigate what could be done with 
state-of-the-art programmable calculators to improve on the 
in situ convergence zone predictions available to air ASW 
units . 

The reader only interested in the calculator programs 
developed may skip immediately to the appendix. 

B. ARE LARGE COMPUTERS NECESSARY? 

In section 5.6 of Ref. 1, Principles of Underwater Sound, 
by R.J. Urick, the author discusses the relative merits of 
two theoretical approaches to obtaining wave equation solu- 
tions in order to describe the distribution of sound energy 
in space and time. Several references are made to the need 
for digital computers to produce sound propagation descrip- 
tions with either theory. Since those comments were made. 
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however, there has been a revolution in the capabilities of 
small programmable calculators. Although it is probably true 
that computers are required to produce a complete description 
of sound propagation in the ocean with one program, calcula- 
tors are capable of solving the different modes of propaga- 
tion one at a time with separate programs to obtain a composite 
description. Examples of simple but fairly adequate calculator 
programs for surface duct, bottom bounce, reliable acoustic 
path, and deep sound channel propagation modes are contained 
in Refs. 2 and 3. These references also contain simple models 
for CZ propagation, but they are based on a mix of ray theory, 
rule-of-thxomb , and empirical data. It was felt that a better 
CZ model needed to be developed. 

C. DESIRABLE CHARACTERISITCS OF A CALCULATOR PROGRAM FOR 

CZ PREDICTIONS 

1. The program should require a minimum of easily avail- 
able input data. The only information not currently available 
but which would be needed by an airborne ASW unit is an SVP 
from the permanent thermocline to the ocean bottom in the 
unit's search area. A chart of this data could easily be in- 
cluded in the environmental package carried aboard the aircraft. 

2. The program should be easy to operate and not require 
the operator to have a great deal of insight into the mathe- 
matical model or the internal operations of the calculator. 

3 . The output should provide ranges to the inner and 
outer edges of all CZ annuli of interest. It should also pre- 
sent the expected TL for all frequencies of interest in each 
annulus . 
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4. The run time of the program should be relatively 
short. This characteristic recognizes that time is impor- 
tant to an on-station unit. 

5. The program should be based on generally accepted 
acoustic theory. This characteristic is desirable because a 
user would probably have more confidence in such a model 
than one based on empirical data and thus applicable only to 

a specific ocean basin. With empirical models, the user often 
wonders if the area he intends to search corresponds to the 
mean set of conditions used to generate the model or is some- 
how different. 

6. Ideally, the program's performance should agree 
closely with the generally accepted large computer models 
currently in use. 
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II. USING I CAPS AS THE STANDARD FOR COMPARISON 



A. RATIONALE FOR ONLY ONE LOCATION PER OCEAN BASIN 

As will be demonstrated, the CZ characteristics of the 
three locations studied vary considerably. The deep sound 
channels which produce CZ phenomena are quite different as 
are the ranges from source to CZ annuli. The objective of 
this study was to produce a mathematical model for CZ pre- 
dictions for use on small programmable calculators. It was 
reasoned that if the model would work for the different con- 
ditions of the three locations studied, it would work for all 
of the variations to be expected within any one of the ocean 
areas . 

B. DESCRIPTION OF ICAPS 

The Integrated Carrier ASW Prediction System (ICAPS) is 
a passive and active acoustic prediction system, developed 
for installation aboard aircraft carriers and other large 
naval vessels which have digital computers. It contains 
four sets of historical environmental data, one each for the 
North Pacific, North Atlantic, and Indian Oceans, and one 
for the Mediterranean and Black Seas. It also contains sev- 
eral production programs for predicting naval sonar system 
performance. The FACT model is used in the passive sensor 
predictions. This is the same model used at FNWC for ASRAPS. 
Reference 4 contains a description of the installation and 
operation of ICAPS in the IBM 360 Computer Center at the 
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Naval Postgraduate School, Reference 5 contains a descrip- 
tion of the mathematics used in the FACT model. 

C. DESCRIPTION OF HISTORICAL ENVI R0NI»1ENTAL DATA FILES 

Figures 1(a) through ICl) depict twelve sound velocity 
profiles produced by ICAPS from its historical environmental 
data files. Figures 1(a) through 1(d) show SVP information 
for the months of February, May, August, and November for 
40N 140W in the Pacific Ocean. Figures 1(e) through 1(h) 
are the same information for 31N 69W in the Atlantic Ocean, 
and likewise. Figures l(i) through ICl) are for 36N 18E in 
the Mediterranean Sea. 

The historical data files used to produce these profiles 
consist of temperature and salinity values for over thirty 
depths, for four seasons of the year, and for many locations 
spaced at one to five degree latitude and longitude inter- 
vals in each ocean area covered. When specific latitude, 
longitude, and date are specified, interpolations are per- 
formed to produce the approximate temperature and salinity 
profiles to be expected at that location and date. This in- 
formation is then converted to an SVP. The output from this 
portion of the system consists of seven columns of values, 
one each for depth in meters and feet, temperature in Cel- 
sius and Fahrenheit, salinity, and sound velocity in meters 
per second and feet per second. The depths associated with 
these quantities begin at ten meter intervals near the ocean 
surface and gradually increase through 25, 50, 100, 250, 500, 
and 1,000 meter intervals as depth increases. The last line 
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of values is for the ocean bottom depth which was part of 
the input data . 

D. COMPARISON OF DEEP SOUND CHANNEL CHARACTERISTICS 

Figure 2 depicts the deep sound channel (DSC) portion 
of the May SVP for the Pacific/ Atlantic and Mediterranean 
coordinates mentioned earlier in a composite graph drawn to 
scale. As can be seen in that figure, DSC characteristics 
of the three areas differ considerably. The vertical ex- 
tent of the channels varies from 1100 meters in the Mediter- 
ranean to 4200 meters in the Atlantic. The change in sound 
velocity between sonic layer depth (SLD) and DSC axis (point 
of minimum velocity) varies from 15 m/sec in the Mediterra- 
nean to 38 m/sec in the Atlantic. The depth of the DSC axis 
varies from 100 meters in the Mediterranean to 1300 meters 
in the Atlantic. Pacific Ocean values are between the others 
for all of those characteristics. Sound velocity near the 
surface is much greater in the Atlantic and Mediterranean 
than in the Pacific, and sound velocity near the bottom of 
the three basins (not shown in the figure) is about 6 m/sec 
greater in the Atlantic than in the Pacific and about 52 
m/sec greater in the Mediterranean than in the Pacific at 
equal depths. Also note the subsurface sound channel lo- 
cated about 100 to 500 meters below the surface in the At- 
lantic profile. 
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E. COMPARISON OF CZ CHARACTERISTICS IN THREE OCEANS 
1 . Method Used in Obtaining Data for Comparison 

There were two primary objectives in gathering 
twelve ICAPS runs from each of the three ocean areas. First, 
it was desired to obtain sufficient data to determine which 
CZ characteristics are common to all areas and which charac- 
teristics are peculiar to specific basins. Secondly, it was 
desired to obtain a standard of comparison for any calculator 
program which might be developed. To fulfill the first ob- 
jective, it was decided to keep the input variables the same 
in all areas, varying them one at a time, in order to better 
compare the differences observed in the various runs. For 
each of the three locations, the inputs varied were season 
of the year and source and receiver depth combination. Re- 
ceiver depths of 60 and 300 feet and source depths of 60 and 
400 feet were used. Each of the ICAPS outputs consisted of 
TL profiles for four frequencies out to a range of 250 kyds . 

Originally, it was intended to collect twelve data 
from each profile. These data were to be the range, width, 

CZ gain, and Transmission Loss for each of the first three 
convergence zone annuli. As it turned out, somewhat less 
data was collected and tabulated. There were several rea- 
sons for this. First, the February SVP in the Mediterranean 
contained no sound channel and therefore no convergence zones 
existed. Secondly, all of the third CZ data for the Atlan- 
tic was thrown out on the grounds that it was almost always 
the same and that it was inconsistent with information from 
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the first two CZ annuli in any particular profile. The 
reason for this occurrence is not know. Finally, it was 
impossible to obtain some of the desired data because of 
the smooth way in which the CZ path blended with other com- 
petitive propagation modes. One could not tell what was 
CZ and what was not in those cases. 

2 . CZ Range and Width Analysis 

The transmission loss profiles produced by ICAPS 
are presented in two formats, a table of TL values for each 
kiloyard of range from the source and a graph of the same 
information. Because the TL values are tabulated at kilo- 
yard intervals, it is impossible to be more accurate than 
that interval in determining where a CZ begins and ends . 
Also, it was difficult to be consistent in picking the 
points representing the edges of CZ annuli because of the 
variety of graph shapes, TL levels, and other propagation 
mode interferences. In any event, an attempt was made to 
satisfy one basic criterion in choosing leading and trailing 
edges of the annuli: Do these ranges best represent the 

apparent location of the annulus regardless of the TL levels 
involved? Admittedly, the ranges picked were often based 
on subjective judgement, and it cannot be stated with com- 
plete certainty that only CZ mode propagation contributed 
to the TL peaks judged to be the CZ annuli. 

Table I contains the CZ range and width data that 
could be gleaned from the ICAPS profiles. In the table, 

RCZi is the range to the inner edge of the first, second, or 
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Table I. ICAPS CZ 
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third CZ annulus to the nearest one half nautical mile. CZW 
is the width of the annulus also the nearest one half nauti- 
cal mile. The third column of numbers is the ratio of CZW 
to the range of the outer edge of the annulus (RCZo) , ex- 
pressed as a percentage. 

After carefully studying this data, the following 
conclusions were made concerning CZ propagation: 

(1) The range to the first CZ is approximately 
14 to 18 nm at the Mediterranean location,' 23 to 27 nm at 

the Pacific location, and 33 to 35 nm at the Atlantic location. 

(2) Range to the CZ decreases and annulus 
width increases as source and receiver get deeper in all 
cases . 

(3) The ranges to the second and third annuli 
are approximately whole number multiples of the ranges to the 
inner and outer edges of the first annulus in all cases. 

(4) The range or width of a CZ annulus does 
not appear to have any significant frequency dependence. 

3 . CZ Gain and Transmission Loss Analysis 

Convergence zone gain is defined as the difference 
between the transmission loss expected under conditions of 
spherical propagation and the actual transmission loss ob- 
served. This definition is expressed in Eq. (1) . 

G = 20 log(r) + a(r) - TL (.1) 

In this equation, G is the CZ gain, r is the range to the 
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CZ annulus, a is the attenuation coefficient associated with 
the frequency of interest, and TL is the actual transmission 
loss observed in the CZ annulus for that frequency. All 
terms in Eq. (1) are in decibels (dB) . 

In actual convergence zones, TL (and therefore gain) 
is by no means a constant value. Contributions of several 
possible propagation paths at any one point and the time 
varying nature of sound paths in the ocean cause coherence 
effects to exist. These effects make TL vary in both space 
and time. Coherence effects are more pronounced at lower 
frequencies (longer wavelengths) where the time varying ef- 
fects are small compared to spatially distributed effects. 

In ICAPS, the more predictable coherence conditions are in- 
cluded in the mathematical model. 

Since a single TL value was desired for the envi- 
sioned calculator model, an attempt was made to pick the 
"average" TL in the ICAPS CZ annuli. As with the range es- 
timates, this called for subjective judgement. Figure 3 
shows a typical ICAPS CZ presentation which has coherence 
effects in evidence. The figure suggests how an "average" 

TL was chosen as best representing that annulus. Two levels 
were chosen (labeled high and low in the figure)' which bracket 
the majority of the TL points within the annulus. The ap- 
proximate midpoint between those levels was then picked as 
"the" TL for that CZ . ^ 

As an extra point of interest, the high and low TL 
levels were studied. It was noted that ICAPS predicts TL 
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Transmission Loss (dB) 




Figure 3 * Fstimating 
Annulus from an ICAF3 



Transmission 
TL Profile. 



Loss 



in a GZ 
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variations from about ± lOdB around the "average” level. If 
this is truly representative of CZ coherence effects, an ASW 
unit armed only with an estimate of the "average" TL in a 
certain CZ annulus should expect to see variations of about 
that magnitude around the estimate in hand. 

Using TL levels estimated by the procedure described 
above, and employing Eq. (1) , CZ gain values predicted by 
ICAPS were obtained and tabulated. The values produced are 
contained in Tables II and III. Table II shows all of the 
data from the Pacific location. Table III contains only 300 
Hz data from the Atlantic and Mediterranean locations. (50, 
850, and 1700 Hz data were omitted from Table III because it 
became obvious during data collection that G is not frequency 
dependent . ) 

Again after careful study, the following conclusions 
were drawn concerning CZ gain : 

(1) CZ gain values range from about eight to 
twenty dB in all three areas observed. 

(2) CZ gain is the same value for first, 
second, and third CZ in any given case. 

(3) In general, CZ gain is independent of 
frequency. An exception to this conclusion is that at low 
frequency (below 300Hz) , especially when source or receiver 
or both are above the SLD and/or near the surface, there is 
apparently somewhat less gain than evident for higher fre- 
quencies. This difference is probably due to stronger dif- 
fraction of the longer wavelengths. 
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C4) CZgain seems to be highest when source 



and receiver are at or near the same depth. 



43 



Month 


Rcvr/Tgt 


Freq 




CZ 


Gain (dB) 




SLD(ft) (ft) 


(Hz) 


1st CZ 




2nd CZ 


3rd CZ 






50 


7 




7 


6 




60/60 


300 

850 


14 

14 




13 

13 


12 

13 






1700 


15 




15 


12 






50 


9 




9 


9 


FEB 

246 


60/400 


300 

850 


12 

11 




12 

11 


11 

12 






1700 


11 




11 


12 






50 


12 




11 


10 




300/400 


300 

850 


12 

10 




12 

11 


11 

12 






1700 


12 




10 


11 






50 


16 




17 


16 




60/60 


300 


16 




17 


16 




850 


16 




17 


16 






1700 


16 




17 


17 






50 


10 




10 


12 


MAY 


60/400 


300 


10 




11 


13 


33 


850 


11 




12 


15 






1700 


9 




11 


12 






50 


14 




14 


14 




300/400 


300 


13 




13 


14 




850 


13 




14 


11 






1700 


12 




11 


13 


Table II. 


ICAPS CZ Gain Data 


Observed 


at 


the Pacific Ocean 




Location . 


(Page 1 


of 2) 
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Month 


Rcvr/Tqt 


Freq 




CZ 


Gain (dB) 




SLD(ft) (ft) 


(Hz) 


1st CZ 




2nd CZ 


3rd CZ 






50 


11 




11 


13 




60/60 


300 

850 


14 

15 




13 

18 


12 

20 






1700 


16 




17 


18 






50 


7 




7 


8 


AUG 

0 


60/400 


300 

850 


9 

12 




12 

14 


13 

15 






1700 


12 




13 


15 






50 


14 




14 


14 




300/400 


300 


11 




14 


12 






850 


11 




12 


15 






1700 


12 




13 


11 






50 


8 




7 


6 




60/60 


300 


9 




12 


12 




850 


15 




18 


16 






1700 


14 




15 


15 






50 


9 




9 


8 


NOV 

98 


60/400 


300 

850 


12 

11 




13 

11 


11 

11 






1700 


11 




11 


11 






50 


11 




11 


11 




300/400 


300 


12 




11 


11 




850 


9 




11 


12 






1700 


10 




9 


11 


Table II. 


ICAPS CZ Gain 


Data 


Observed 


at 


the Pacific Ocean 




Location. (Page 2 


of 2) 
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Ocean 



Month 

SLDCft) 



Rcvr/Tgt 

(ft) 



1st CZ 



300 Hz CZ Gain 



2nd CZ 



3rd CZ 



FEB 

328 



60/60 


17 


14 


60/400 


13 


13 


300/400 


17 


18 



u 

H 



Eh 



Eh 

< 



< 



W 



X 

W 

Eh 



M 

Q 

W 



MAY 

n 


60/60 


19 


14 


60/400 


16 


12 


u 


300/400 


13 


17 


AUG 

n 


60/60 


15 


17 


60/400 


10 


13 


U 


300/400 


14 


15 


NOV 

n 


60/90 


13 


11 


60/400 


12 


10 


U 


300/400 


13 


12 



MAY 

10 


60/60 

60/400 

300/400 


14 

12 

17 


13 

10 

18 


16 

10 

20 


AUG 


60/60 


17 


19 


17 


0 


60/400 

300/400 


11 


9 


10 


NOV 


60/60 


11 


13 


15 


XT' 


60/400 


10 


11 


12 


300/400 


12 


14 


15 



Table III. ICAPS 300Hz CZ Gain Data Observed at Atlantic 
and Mediterranean Locations - 
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III. CZ RAY THEORY ANALYSIS AND MODEL DEVELOPMENT 



A. CZ RANGE AND WIDTH 

It was decided to use ray tracing as the method for de- 
termining CZ range and width because of the simplicity of 
the mathematics involved and because of the intuitive appeal 
of sound rays depicting the propagation of sound. The alter- 
native approach, that of normal mode theory, was rejected on 
the grounds that it would be much more complicated, requiring 
capabilities far beyond those available in the calculators 
at hand. ' ’■ 

Figure 4 shows four sound rays of particular interest in 
CZ propagation. The order of these rays is described for the 
"typical" case in the following discussion. An "atypical" 
case will be mentioned later. 

Ray #1 departs the SLD at zero degree depression angle. 

It reaches its greatest depth at the bottom of the DSC and 
returns to the SLD at some particular range and at horizontal 
incidence. The horizontal range from SLD to SLD is termed 
cycle distance. The cycle distance for this ray is desig- 
nated rQ . 

Ray #2 is the next ray of interest found as the departure 
angle from the SLD is increased downward. This ray passes 
down through and below the bottom of the DSC before turning 
back, upward. It returns to the SLD at the shortest range 
from the starting point of any ray within the bundle of rays 
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Sound Velocity |- Range 
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undergoing CZ refraction. This range is designated 
The angle of departure for this ray is designated 
Each ray between #1 and #2 crosses all of the previous (les- 
ser departure angle) rays on its way up from its lowest 
depth. 

As departure angle from the SLD is further increased, 

the next sound ray of interest, #3, is located. This ray 

has a cycle distance equal to r^ . Its maximum depth is 

greater than that for ray #2. It departs from and arrives 

back at the SLD at an angle designated 9 . Rays between 

^ ^ rswp ^ 

#2 and #3 do not cross each other, but they do cross the 
earlier rays on their way back up to the SLD. 

In the CZ annulus, the rays between #1 and #2 sweep in- 
ward toward the source as departure angle increases. After 
ray #2 they sweep out away from the source as angle increases 
further. For this reason, the region formed by rays between 
#1 and #3 is called the reswept zone . 

Finally, as angle of departure from the SLD is increased 
to maximum angle for CZ propagation, we observe ray #4. This 
ray turns upward at a depth equal to the water column depth 
at that location. It returns to the SLD at the greatest dis- 
tance of all CZ refracted rays. Rays departing the SLD at 
angles greater than that for ray #4 would be reflected off 
the bottom and are not of significance for the CZ propagation 
path. 

As mentioned earlier, this progression of rays exists 
in a "typical" CZ situation. If, however, the ocean bottom 
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were more shallow, cycle distance for ray #4 would be re- 
duced. If the bottom were shallow enough, ray #4's cycle 
distance would be less than r^ . In that case, the reswept 
zone would be reduced to the region between rays #2 and #4. 
This situation is called the "atypical" case. 

B. RANGE AND WIDTH MODEL 

In the calculator programs developed, provision is made 
for entering and storing a five point sound velocity profile 
which defines the DSC only. The first depth and velocity 
pair entered (Dj^,Cj^) equate to the appropriate values found 
at the SLD. The fifth depth entered (D^) is the depth at 
the bottom of the DSC where sound velocity is equal to that 
at the SLD. The other three depth/velocity pairs must be 
picked subjectively from a graph of the SVP of interest. If 
a mixed layer exists, the gradient in that layer is taken to 
be 0.Q2 sec ^ (a purely pressure induced gradient to two 
place accuracy) . The program calculates the four layer gra- 
dients within the DSC profile entered, and uses the fourth 
(deepest) layer gradient in ray calculations that occur below 
the DSC. It would have been desirable to allow several more 
points in the SVP, but calculator data storage capacity and 
program step limitations preclude more than five depth/veloc- 
ity pairs. 

The overall scheme used to predict CZ range and width is 
to trace a series of rays starting at the SLD with a zero de- 
pression angle ray. That first ray yields rQ which is stored. 
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Then an iterative process is begun in which the angle is in- 
cremented and each succeeding cycle distance determined is 

compared to the previous one until r . and 9 . are found. 

min rmin 

C6rmin stored for use in the CZ gain and TL program, to 

be discussed later.) Corrections are then made to r . and 

min 

rQ to account for surface duct effects (if any) and source 
and receiver depth separation from the SLD. Range to the 
inner edge of the CZ is plus corrections, and range to 

the outer edge is r^ plus corrections. 

In the first attempt to produce a calculator program, 
the cycle distance of ray #4 (the ray just grazing the bot- 
tom) was compared to r^ . The greater of the two was picked 
as the basic distance for determining range to the outer 
edge of a CZ . Later on, this portion of the program had to 
be deleted to save program steps. The final programs de- 
veloped ignore bottom depth and do not include rays outside 
the reswept zone in determining annular width. This is prob- 
ably a shortcoming of the programs but the seriousness of 
the errors it causes will not be known without further study. 

A commonly applied rule-of-thumb states there must be 
a minimum 300 fathoms of depth excess (water column below 
the DSC) in order to have "reliable" CZ conditions. It was 
observed that a fully developed reswept zone existed in 
every case in the locations studied, and separate calcula- 
tions showed that somewhat less than 300 fathoms depth ex- 
cess was required to complete the zone. Therefore, a program 
user should consider the 300 fathom rule-of-thumb before 
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running the range and width program. With less than 300 
fathoms depth excess, the possibility exists for an "atypical" 
CZ propagation situation where the reswept zone is reduced in 
width due to bottom ray limiting. 

Another program shortcoming involves an assumption that 
both source and receiver would be more shallow than the sec- 
ond DSC SVP point chosen (depth - In other words, the 
programs were designed to allow for source/receiver depths 
within the mixed layer or the first isogradient layer below 
the SLD. After the five point SVP is entered and the grad- 
ients computed, source and receiver depths are entered and 
converted to velocities. The programs determine these veloc- 
ities (C_ and C ) by subtracting an appropriate amount from 
the velocity at the SLD, The amount subtracted is determined 
by depth separation from the SLD and by the gradient in either 
the ML or the first layer below the SLD. If source or re- 
ceiver depth is greater than D^t sound velocity should be de- 
termined by correcting C2 (the velocity at D2) and by using 
g2 (the second layer gradient) , Since this is not done, ve- 
locities for source/receiver depths below D2 will be in error 
(usually too low) . Source and receiver velocity errors are 
carried over into Ar„ and Ar range corrections. If the 
velocities are too low, the range corrections will be too 
large. This is normally a rather insignificant source of 
total range error, however, since Ar^ and Ar^^ errors will be 
a small fraction of the magnitude of those terms and because 
the range correction terms are small to begin with. 



The mathematics of ray tracing in isogradient layers is 
quite straightforward. By Snell's Law, 

C C 

1 _ 2 _ (A constant . 

cos 9^ cos 02 for each ray) ^ 

the angle of a ray departing a layer can be determined from 
the angle of entry into that layer. In Eq. (2) , is the 
sound velocity where the sound ray enters a layer, 9^ is the 
angle of entry, C 2 is the sound velocity where the ray de- 
parts the layer, and 02 is the angle of departure. 

Rays travel in circular arcs within constant gradient 
layers, and the radius of curvature is: 



R 





(3) 



and 9^ are as defined above and g^ is the gradient within 
the layer (in this case, layer 1) . 

Finally, the horizontal distance traveled by a ray 
while traversing a layer is : 



Ar = 



R (sin 9^ - sin 9^) 



( 4 : 



Figure 5 demonstrates an example application of Eqs . 

2 through 4. It should be noted that absolute value signs 
are used in Eqs. 3 and 4 because the gradient in Eq. 3 and 
the difference of sines in Eq. 4 may be positive or negative, 
while R and Ar are always positive. 

The programs use these equations to compute the horizon- 
tal range increments each ray accumulates within the four 
layers, doubles each term Cto account for the downward and 
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Depth 




Figure 5* Horizontal Distance Traveled Within 

an Isogradient Layer. 
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upward passes through each layer) , and then sums the terms 
to obtain cycle distances. The fourth layer requires a 
slightly different treatment because the rays become horizon- 
tal and then turn back upward within that layer. Essentially 
the same formulas are used, however. The equations are also 
used to compute the range correction terms. 

In considering the various possible ray paths between 
source and receiver, it was decided there were four basic 
situations which could occur: 

1) No mixed layer, both source and receiver 
below the SLD. (.Deep/Deep) 

2) Mixed layer present, both source and receiver 

below the SLD. ( Deep/Deep/ML) 

3) Mixed layer present, both source and receiver 
within the layer. (Shal/Shal) 

4) Mixed layer present, source or receiver above 
the SLD, the other below. (Crosslayer) 

The only difference between the first two cases is the mixed 

layer effect in case 2. That effect causes a widening of 

annuli due to spreading of sound rays as they travel up to 

the surface and back down to the SLD within the layer. The 

mixed layer effect is also included in the third and fourth 

cases above. 

It was originally intended to include all four cases in 
one range prediction program. Again due to calculator limi- 
tations, it was necessary to use two programs to cover the 
four possibilities. The first range program (labeled Deep/ 
Deep) is for cases 1) and 2) above when both source and re- 
ceiver are below the SLD whether or not an ML exists. The 
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second range program (labeled Shal/Shal or Crosslayer) is 
for use in cases 3) and 4) above when source or receiver or 
both are above the SLD. 



Formulas used to determine range to inner edge of the 
first CZ CRCZi) and range to the outer edge of the first CZ 
(RCZo) follow: 



RCZi = 


r . 
min 


- ^>^3 - 






Deep/Deep 


= 


r . 
min 


+ Arg + 






Shal/Shal 


= 


r . 
min 


+ Ar 

s 


^^R 




Crosslayer 


RCZo = 






+ Arg 


+ Ar^ 


Deep/Deep 

Deep/Deep/ML 


= 




2 ArQ + 


- 




Shal/Shal 


= 


^0 + 


2 Ar^ + 


Arg + 




Cross layer 


these equations , r . 

^ min 


and rQ 


have been 


previously def 



Ar and Ar are the respective horizontal range corrections 
which account for source and receiver depth separation from 
the SLD, and 2 Ar^ is the correction for mixed layer effect. 
Figures 6(a) through 6(f) (not to scale) depict the RCZi and 
RCZo formulas in graphic form. Ranges to second and subse- 
quent CZ annuli are taken to be integer multiples of the 
ranges produced. 

Two items of interest, both evident in Figs. 6(a) - 6(f) , 
are worth mentioning at this point. First, acoustical reci- 
procity is envoked and the more shallow of source and re- 
ceiver is always treated as "source" of the sound rays within 
the calculator programs. Secondly, only those sound rays 
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Figure 6(a). RCZ^ for Deep/Deep Case. 
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which experience no more than one ocean surface reflection 
between source and receiver are considered in this model. 

Both of these conventions are commonly applied to ray trac- 
ing models. Although they theoretically have' little or no 
effect on model results, they greatly simplify the work of 
programming a ray tracing model. 

C. CZ GAIN AND TRANSMISSION LOSS MODEL 

In general, transmission loss is defined as ten times 
the logarithm of the ratio of sound intensities measured at 
one meter from a source and at range r from that source. 

TL = 10 log 

Intensity has units of power per unit area. The change in 
intensity between one meter and range r is due to geometric 
spreading of the power over a different amount of area and 
due to attenuation of some of the power through absorption, 
scattering, diffusion, etc. 

In ray tracing theory, it is assiimed there is no sound 
power transfer across sound rays. Therefore, the power flow- 
ing from a source between two sound rays remains between 
those rays and travels out in a direction parallel to the 
ray paths. Determining the portion of transmission loss due 
to geometric spreading (TL^) under this assumption reduces to 
finding ten times the logarithm of the ratio of areas (at 
range r and at one meter) penetrated by the power between 
the two rays perpendicular to the direction of travel. 
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TLg = 10 log ^ 

A mathematicai development of this technique is contained on 
pages 119-121 of Ref. 1. 

Figure 7 shows how this method was adapted for use in 
the CZ Gain and Transmission Loss portion of the calculator 
model developed. The area CA^) at one meter from the source 
is the product of area height and area circumference. The 
sound rays bounding the area above and below are the minimum 
and maximum departure angle rays which produce the reswept 
zone in the CZ annulus. The angular spread of those rays 
CA0) in radian units times the sphere radius (1 meter) is 
the area height. Cosine of the average angle of departure 
of the rays (0^^) times the sphere radius times 2tt is circum- 
ference of the area. Therefore: 

A^ = 27t A9 cos 0^ (m^) 

In the CZ, the area (^2^ over which the same power is 
distributed is also found by a product of area circumference 
and area width. Circumference is 2'n' times range to the CZ 
(RCZi) . Width of the area perpendicular to the sound rays 
is the product of CZ annulus width (CZW) and the sine of the 
average angle of arrival of the rays at the receiver depth 
(02) . Therefore; 

2 

A 2 = 2ir RCZi CZW sin 0 2 Cm ) 
and the geometric TL expression becomes : 
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I'igure 7. Geometric Transmission Loss by Kay Tracing Method 



RCZi CZW sin 0_ 

TL = 10 log T -5 s 

g ^ AS cos 

Substituting this expression back into Eq. 1, which is 
the definition of CZ gain, and reducing to simplest form 
yields the algorithm used to determine G in the calculator 
model : 

RCZi A0 cos 0, 

^ cz - W - sTh"^ ^ ^ 

To implement this algorithm, the program has only to 
determine the angular terms since RCZi and CZW are available 
from the range program results. After one of the range pro- 
grams has been run, the user loads the G and TL program into 
calculator memory without altering the contents of the data 
storage registers left from the range program. Then the 
iterative ray tracing process begun in the range program is 
continued in the gain program. The angle of departure of 
sound rays from the SLD is incremented beyond 0 (left in 

storage from the range program) and cycle distances produced 
are checked for approximate equality with r^ . In this way, 
the ray which completes the reswept zone is found, and its 
angle of departure from the SLD is 0 . Then 0 and 

zero degrees (the angle for the ray producing cycle distance 
r^) are converted to angles of departure from the source 
depth, 0„ and 0^^ respectively, and angles of arrival at 

the receiver depth, 0 and 0 respectively, using Snell's 

RR RO 

law. The angular terms in Eq. 5 are then computed using 
the following formulas: 
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II 

CD 

cn 


+ 


CD 

cn 

0 
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- 


®so 
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- 


0 

O) 

CD 
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- ®SR 


+ 


CD 

cn 

0 


Shal 
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2 



fl _ ®R0 ®RR 
®2 2 

Recall that "source" in the model refers to the more shallow 
of source and receiver. Therefore, the deep "source" forms 
of these formulas are used only after using the Deep/Deep 
range program. In all other cases, the "source" is consid- 
ered to be shallow. Figure 8 depicts these angular relation- 
ships for the various depth conditions. 

It should be pointed out there are two inherent errors 
in the angular quantities determined. First, the possible 
source and receiver sound velocity errors mentioned earlier 
could cause the gain algorithm angles to be slightly off. 

This would only occur if depths greater than D 2 were entered 
for source or receiver or both. Secondly, 9 is found 

for the ray which has cycle distance equal to r^ at the SLD. 
Since the actual CZ ray bundle departs from the source depth 
(vice SLD) and arrives at the receiver depth (vice SLD) , the 
ray which completes the reswept zone will probably be differ- 
ent than the ray used and it will have a slightly different 
angle crossing the SLD. These angular errors will cause the 



63 



0 ° 






r’igure 8. Determining angular terms for 

CZ Gain algorithm. 
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greatest CZ gain error in the sin 62 term of the algorithm. 
Since sine is directly proportional to angle at small angles, 
an error of a factor of two in quite possible event) 

could cause a gain error of approximately 3 dB . 

Once CZ gain is computed and stored, the sound frequency 
of interest is entered, and the attenuation coefficient is 
calculated using Thorpe's equation (p. 102, Ref. 1): 

(dB/m) ( 6 ) 

In Eq. 6 , f is in kHz, and the constant in front of the ex- 
pression converts attenuation coefficient from dB/kyd to 
dB/m. The program user enters frequency in Hz, and the pro- 
gram performs the conversion to kHz. 

til 

Finally, the transmission loss in the n CZ annulus 
for the frequency of interest (TL^) is determined by: 

TL^ = 20 log (n RCZi) + z (n RCZi) - G (7) 

til. 

In this equation, the subscript n denotes the n CZ annulus, 
the range to which is n times RCZi. 

After a range program is run, and after the gain portion 
of the G and TL^ program has been completed, TL^ values for 
a variety of frequencies and CZ annuli may be rapidly ob- 
tained for the SVP , source depth, and receiver depth condi- 
tions entered. If, however, a different set of source/receiver 
depth conditions are also of interest, the entire procedure 
beginning with the appropriate range program must be per- 
formed again. 



a = (0.001094) 



0.1 r 

1 + f“ 



40 f^ 
4100 + f 
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D. CALCULATOR PREDICTIONS COMPARED TO ICAPS 



1 . Choosing the SVP Points for the Program 

The five point SVP limitation of the ray tracing 
procedure is a rather serious handicap in many situations. 
Actual sound velocity profiles not only are curvilinear in 
overall shape but also have many small scale features and 
they are time varying functions . Approximating these curves 
with only four straight line segments presents a difficult 
challenge . 

In general, matching the gradients, sound velocities, 
and associated depths are all important in choosing SVP 
points. The greatest potential for causing large range pre- 
diction errors occurs when the SVP contains an extensive 
near surface layer with very slight velocity gradient. The 
horizontal distance traveled by a shallow depression angle 
ray within such a layer varies considerably with small 
changes in the gradient or layer thickness. Under such con- 
ditions, then, it is extremely important to match those 
characteristics as closely as possible. 

Another important item to carefully match is the 
sound velocity at the DSC axis. This velocity determines 
the maximum angle of depression for each ray prior to com- 
mencing upward refraction. The horizontal distance traveled 
by a ray below the axis is highly dependent on that angle. 

Table IV contains the five point sound velocity 
profiles picked by the author for use in comparing the pro- 
gram performance to ICAPS predictions. Depths in the table 
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5 Points 


FEB 


MAY 


AUG 


NOV 
















PACIFIC 






D1 


75 


10 


0 


30 


Cl 


1489 . 8 


1500.2 


1508.6 


1497 . 7 


D2 


340 


80 


100 


100 


C2 


1480.0 


1490.0 


1483.0 


1484.0 


D3 


500 


600 


600 


600 


C3 


1476.2 


1477.0 


1475.5 


1476.0 


D4 


1060 


1600 


1750 


1400 


C4 


1480.0 


1485.0 


1486.5 


1483.0 


D5 


1940 


2610 


3120 


2460 






ATLANTIC 






D1 


100 


0 


0 


0 


Cl 


1524.6 


1530.9 


1541.1 


1535.5 


D2 


550 


100 


125 


125 


C2 


1522.5 


1522.0 


1522.5 


1522.5 


D3 


1080 


575 


550 


625 


C3 


1493.0 


1524.0 


1524.5 


1524.5 


D4 


1750 


1225 


1200 


1175 


C4 


1495.0 


1488.5 


1487 . 5 


1489.0 


D. 

D 


3780 


4175 


4760 


4460 




MEl 


DITERRAMEAN 






D1 




10 


0 


10 


Cl 




1526 . 6 


1537.6 


1528.7 


D2 




30 


50 


50 


C2 




1518.0 


1517 . 5 


1517.5 


D3 




100 


100 


125 


C3 




1513 .0 


1513.0 


1512.0 


D4 




800 


700 


900 


C4 




1521.5 


1519 . 5 


1522.5 


D5 




1125 


1800 


1290 



Table IV. Five Point Sound Velocity Profiles. 

(Depths in meters, velocities in m/sec) 
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are in meters, and sound velocities are in meters per second. 
The reader may want to plot these points on the graphs of 
Figs. 1(a) through 1(1) so he may see how the four isogradi- 
ent layers picked match the ICAFS profiles. It should be 
pointed out that only the initial selection of SVP points 
was used in the subsequent comparisons of calculator model 
results to ICAPS predictions. Since an ASW aircrewman using 
the programs in attempting an in situ prediction of acoustic 
conditions would not be able to judge whether SVP point ad- 
justments would improve or degrade prediction accuracy, it 
was felt that comparing resuts of the initial SVP point se- 
lection with ICAPS would be more meaningful to the objective 
of developing the calculator programs. 

Comparing calculator model predictions to ICAPS pre- 
dictions in a definitive statistical manner was not done. 

The main reason for this was alluded to in the preceding para- 
graphs. Since the SVP points entered in the calculator pro- 
gram must be picked subjectively by the person using the 
program and since it is unlikely different people would pick 
the exact same points off any given SVP, it is clear that 
calculator results can be expected to vary from operator to 
operator . 

2 . CZ Range and Annulus Width Comparisons 

The calculator range programs produce one value each 
for RCZi and RCZo for any given SVP, source depth, and re- 
ceiver depth situation. Under the same set of conditions, 
ICAPS yields four sets of RCZi and CZW values, one set for 
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each of the four frequencies entered. In order to compare 
the calculator performance to ICAPS it was first necessary 
to reduce the ICAPS predictions to one value each for RCZi 
and CZW for each SVP/source/receiver condition. This was 
done by simple averaging to eliminate the frequency variable 
from the ICAPS range and width predictions. 

Figures 9(a) - 9(c) display range and width compari- 
sons in graphical form for the Pacific, Atlantic, and Medi- 
terranean locations respectively. In each figure the double 
barred lines represent the ICAPS first CZ annuli predictions 
(averaged over frequency) , and the single barred lines repre- 
sent the calculator predictions. Numerical values for inner 
and outer first CZ ranges may be obtained from the scales at 
the tops of the figures. 

In all, there were 32 cases where these graphical 
comparisons could be made. The following comments pertain 
to those comparisons: 

a) In 30 of the 32 cases the calculator annuli 
overlap at least a portion of the ICAPS annuli. 

b) In 14 of the 32 cases the calculator annuli 
are completely contained within the limits of the ICAPS 
annuli . 

c) In all 32 cases RCZi ranges predicted by 
the calculator were greater than those predicted by ICAPS . 

In the 12 Pacific cases, the calculator RCZi values were ap- 
proximately 1.7 nm greater than ICAPS on the average. In 

the 12 Atlantic cases, the average difference was approximately 
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Figure 9(a.) . Comparisons of ICAPS and Calculator CZ 
Annuli Predictions. Pacific Ocean. 
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Tgt/Rcvr 

(ft) 35 ^0 Range (nm) 




Legend: ICAFS < i Calculator program i 1 

Figure 9(b) . Comparisons of ICAFS and Calculator CZ 
Annuli Predictions. Atlantic Ocean. 
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Tgt/Rcvr 
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60/60 
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Range (nm) 
20 



I 60/400 
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300/400 
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I 1 



I 60/hoo 

300/400 



60/ 60 



0) 

s 6o/4oo 

<u 



300/400 



Legend; 
Figure 9 (c) . 



ICAPS " I Calculator program 1 1 

Comparisons of ICAPS and Calculator CZ 
Annuli Predictions. Mediterranean Sea. 
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2.8 nm. And in the eight Mediterranean cases, 1.5 nm was 
the mean difference. 

d) In 27 of the 32 cases the CZ width predic- 
tions from the calculator were more narrow than the ICAPS 
predicted widths. Three of the five cases where calculator 
CZW exceeded ICAPS CZW were from the February SVP in the At- 
lantic location. That SVP contained a very deep (500 meter) , 
nearly isovelocity layer near the surface. In such a pro- 
file, CZ refraction produces ray paths that are spread over 
a very wide (in this case 16-20 nm) annulus. Only the rays 
which return to the SLD within the first few nm at the inner 
edge of that annulus experience sufficient convergence to 
produce detectable CZ gain, however. Going from inner to 
outer edge of such an annulus the CZ refracted rays rapidly 
fan out experiencing progressively less convergence and pro- 
ducing progressively less CZ gain. Additionally, if the 
bottom grazing ray were considered, it would be seen to limit 
the reswept region of this type annulus to something far less 
than that indicated. Since the calculator model fails to ac- 
count for either of these factors, it fails rather dramat- 
ically to produce a "practical" CZ annular width from this 
SVP type. 

In summary, the calculator model produces CZ annuli 
that roughly agree with those produced by ICAPS in all three 
ocean basins considered. Calculator RCZi ranges are 5-10% 
greater on the average than the ICAPS values . Calculator 
CZW predictions (excluding the Atlantic February SVP) are 



73 



40-50% narrower than ICAPS widths on the average. And, the 
Atlantic February case indicates there is at least one SVP 
type in which the calculator model fails to produce even 
marginally acceptable results for CZW. 

3 . CZ Gain Comparisons 

As with CZ range and width comparisons, it was nec- 
essary to average the ICAPS gain data with respect to fre- 
quency before calculator gain predictions could be compared. 
The estimated ICAPS gain values in Tables II and III were 
thus reduced to one number for each SVP, source, and receiver 
condition. Table V, CZ Gain Prediction Comparisons, contains 
numbers that represent the difference between calculator gain 
predictions and the averaged ICAPS values. Minus signs in 
the table indicate those cases where calculator gain was less 
than the ICAPS value. 

As with the range and width comparisons, the worst 
agreement occurred in the Atlantic winter SVP case. Since 
CZW is a term in the gain algorithm, the extremely wide an- 
nuli predicted by the calculator caused gain values to be 
far too low for the three source/receiver conditions assoc- 
iated with that SVP . 

Excluding the Atlantic winter SVP case, the follow- 
ing comments can be made concerning the other 29 CZ gain 
comparisons : 

a) Calculator gain values ranged from 7.3 dB 
lower to 5 dB higher than the averaged ICAPS values. 
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^ICAPS ^ ^^calc 



^ICAPS^ 



Source/ 



SVP 

Profile 


Rec eiver 
(ft) 


Pacific 


Atlantic 


Mediterranean 




60/60 


13.4/-5.4 


15.5/-IO.5 




FEB 


6 o/ 4 oo 


11 . 4 /- 2 .4 


13.0/-7.0 






300/400 


11 . 2 /- 2. 2 


17.5/-II.5 






60/60 


16 . 4 / 0 . 6 


I6.5/-O.5 


14 . 3 / 3.7 


MAY 


6 o/ 4 oo 


11 . 6 /- 2. 6 


14 . 0 /- 2.0 


IO.7/-I .7 




300/400 


12 . 7/-0. 7 


15.O/-3.O 


18.3/-6.3 




60/60 


15.8/ 2.2 


16.0/ 2.0 


17.7/-O.7 


AUG 


6 o/ 4 oo 


12 . 8 /- 1 . 8 


11.5/ 0.5 


10.0/ 0.0 




300/400 


12 . 4 / 0.6 


14 . 5 /- 1 .5 






60/60 


l 4 . 0 /- 4. 0 


12.0/ 5.0 


13.0/ 4.0 


NOV 


60/400 


II.3/-7.3 


11.0/ 1.0 


11 . 0/-3. 0 




300/400 


10.7/ 1.3 


12.5/ 0.5 


13.7/-2.7 



Table V. CZ Gain Prediction Comparisons. 
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b) In 22 of the 29 comparisons, calculator 
values were within 3 dB of ICAPS. 

c) In nine of the 29 comparisons calculator 
values were within one dB of ICAPS. 

d) On the average, calculator gain values were 
approximately one dB less than ICAPS. This result is incon- 
sistent with calculator CZW results in light of the gain mod- 
el used. Since calculator CZW values averaged only slightly 
more than half the ICAPS widths, it would have been more con- 
sistent if calculator gain values turned out two to three dB 
higher than ICAPS (acoustic power being spread over less area 
in the CZ annuli, other things being equal) . Perhaps an ex- 
planation for this apparent discrepancy is that the calcula- 
tor model does not consider the contribution of surface 
reflected energy adding to the energy from upward traveling 
sound rays at the receiver depth. In an actual CZ annulus 
the downward traveling, surface reflected energy adds approx- 
imately three dB to the CZ gain over much of the annulus 
width. Apparently, the FACT model in the ICAPS system in- 
cludes this consideration. It is also apparent that neglect- 
ing surface reflected energy in the calculator gain model has 
the effect of canceling errors that should result from CZW 
values being too narrow. 

In summary of the gain results , it can be said that 
the ray tracing technique used in the calculator model worked 
reasonably well. Since three fourths of the comparison cases 
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resulted in gain values within three dB of the estimated 
ICAPS figures, TL values from the calculator displayed the 
same close agreement. 
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IV. CONCLUSIONS 



A. LIMITATIONS OF THE MODEL DEVELOPED 

The HP-67/97 calculators used in programming the CZ pre- 
diction model were stretched to their limits in both data 
storage and program step capacity. Although not known for 
certain, the author feels significantly more accurate results 
would be possible from a calculator with only moderately 
larger storage capacity. 

The data storage limitation which allowed only five SVP 
points to be entered is quite restrictive and no doubt plays 
a large role in the CZ range and width inaccuracies obtained. 

Program step capacity forced several short cuts to be 
taken which again would not have been necessary with a mod- 
erately larger program memory. Two separate range and width 
programs were required due to insufficient program space to 
incorporate tests for different source and receiver depth 
cases. Also, source and receiver depths are strictly allowed 
only within the upper two SVP isogradient layers because pro- 
gram space was not available to check for the correct layer 
if all depths were allowed. Additionally, and perhaps the 
greatest source of CZW errors observed, program step limita- 
tion prevented incorporating a method of considering the 
bottom limited CZ sound ray in determining the range to the 
outer edge of a CZ annulus. The program developed ignores 
the bottom entirely and considers only the reswept zone in 
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predicting annular width. Since calculator CZW results were 
considerably shorter than those indicated by ICAPS, it is 
assumed the discrepancy is due to not considering CZ rays 
beyond the reswept zone. The first priority in making im- 
provements to the calculator model, should a larger capacity 
machine be implemented, would be incorporating a better meth- 
od for selecting the ray which defines the outer limit of 
the CZ annulus. 

B. USEFULNESS OF THE MODEL DEVELOPED 

The degree of success in producing a useable CZ prediction 
model for handheld calculators must be determined by consid- 
ering the objectives set forth in the first section of this 
study. The central idea was to ascertain if a calculator 
model would improve on ASRAPS CZ prediction accuracy in the 
case where BT conditions determined jui situ differed from 
those used to generate the ASRAPS TL profiles. Inherent in 
this objective is the assumption that ASRAPS TL profiles 
generated primarily from climatological data would be in 
error due to lack of input data accuracy. Also inherently 
assumed is that given identical input data the calculator 
model would produce less accurate results than the digital 
computer model(FACT) used in ASRAPS(and ICAPS) due to obvious 
differences in data and program capacities. The real question 
then is a trade-off comparison: Will the basically less accu- 

rate calculator model produce better CZ predictions with actual 
environmental data than the more sophisticated digital com- 
puter model which had only climatological input data? 
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Before addressing the answer to this question, character- 
istics of the calculator model developed will be compared to 
the list of six desirable characteristics described in sec- 
tion I . 

1. Easily available input data. 

The data required are an SVP, assumed source 
depth, hydrophone depth, and frequency of interest. The 
only portion of this information not presently available to 
ASW aircrews is that part of the SVP below the 1,000 ft 
depth limit of the AN/SSQ-36 bathythermograph buoy. SVP data 
from the surface to 1,000 ft (the area where seasonal and 
diurnal variations predominantly occur) is easily obtained 
from the BT buoy information. 

2. Ease of program operation. 

Anyone familiar with HP-67/97 calculator use 
could operate this program without additional training. 

3. Output data. 

The program provides CZ annulus range and width 
as well as TL values for all frequencies of interest in all 
annuli of interest. 

4 . Short run time . 

To run a complete program requires approximately 
10 minutes once SVP data is obtained. Deploying a BT buoy 
and converting the temperature trace to an SVP would take an 
additional 10-15 minutes. 
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5. Based entirely on acoustic theory. 

The program uses only ray tracing techniques in 
producing its output terms. 

6. Agreement with large computer models. 

Calculator CZ ranges obtained averaged 5-10% 

greater than ranges obtained from ICAPS . CZ width results 
averaged only 40-50% of those obtained from ICAPS. And, 
there was one SVP case studied (winter, Atlantic) in which 
the calculator CZW results were very different from ICAPS. 
That SVP case was considered a failure of the calculator 
model, and it must be conceded the model does not work for 
all CZ situations. Excluding the obvious CZW failure SVP 
case, TL values from the calculator averaged about one dB 
lower than ICAPS with extreme deviations observed ranging 
from -7.3dB to +5dB around the ICAPS values. Additionally, 
calculator results can be expected to vary from operator to 
operator since SVP points must be picked subjectively from 
an SVP graph . 

Returning to the main objective of the study, the author 
feels that only half of the trade-off question has been 
answered. An easily operated, purely theoretical model was 
developed which works for most (but not all) CZ producing 
SVP conditions. And a measure of its accuracy compared to 
the sophisticated FACT computer model was obtained. Yet to 
be answered is how inaccurate ASRAPS CZ predictions are when 
observed BT conditions differ from climatological conditions. 
This portion of the question is very difficult to answer and 
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in fact would be a very large study in itself. Generally, 
inaccuracies must range over a scale from insignificant to 
considerable as environmental deviations range from slight 
to great. The most likely variables affecting degree of in 
accuracy are surface water temperature, mixed layer depth, 
and in layer and below layer gradients. The effects of vary 
ing these or other possible factors one at a time or in 
various combinations on CZ range, width, and gain must be 
known before the entire question can be answered. Further, 
definite magnitudes of environmental factor deviation must 
be determined so that a person can judge when ASRAPS inaccu 
racies are likely to be greater than the calculator model 
inaccuracies. Until these points are answered it would be 
inappropriate to recommend use of the calculator model as a 
routine method of updating ASRAPS CZ predictions in situ. 
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APPENDIX 



HP-67/97 Calculator Programs for Convergence 
Zone Range, Width and Transmission Loss Predictions 



Steps required to use the programs : 

1. Deploy a bathythermography buoy in the operating area 
of interest. 

2. Convert the BT buoy information to a sound velocity 
profile of the upper 1,000 ft of the ocean area. 

3. Combine the upper SVP data with a graph of climato- 
logical SVP data which depicts sound velocity conditions 
below the 1,000 ft level. 

4. Pick five points from the combined SVP graph which 
best represent the deep sound channel portion of the SVP. 

The first point should be at the sonic layer depth, the fifth 
point at the bottom of the DSC where sound velocity equals 
that at the SLD, and the other three points at points on the^ 
graph such that when straight lines are drawn to connect the 
five points they create a linearly segmented SVP which matches 
the actual SVP as closely as possible. 

5 . Pick the appropriate CZ Range and Width program to be 

used as follows: If both source and receiver are below the 

SLD use the Deep/Deep program. If source or receiver or both 
are above the SLD use the Shal/Shal or Crosslayer program. 

6. Load and run the appropriate Range and Width program 
according to the accompanying instructions. 

7. Leaving the calculator power on and data storage reg- 
isters unchanged, load and run the Gain and Transmission Loss 
program according to its instructions. 



A word about units : 

As currently written, the programs use metric units; 
meters for depths, and m/sec for sound velocities. To con- 
vert the programs for english unit input data, feet for 
depths and ft/sec for velocities, the conversion factor 
1,852 m/nm should be changed to 6,075 ft/nm where occurring. 
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User Instructions for CZ Range Programs 



Step 



Instructions 



Input Keys Output 



1 



2 

3 



4 



Enter 5-point SVP which 
describes the Deep Sound 
Channel. (D^=SLD, D^= 

bottom of DSC.) 




A 

A 

A 

A 

A 

A 

A 

A 

A 



D 

C 

D 

C 

D 

C 

D 

C 

D 



1 

1 

2 

2 

3 

3 

4 

4 

5 



Press R/S to perform pre- none R/S 

liminary calculations 



Enter Receiver Depth and D^^ 

Source Depth. 

^S 



ENTER 

R/S 



Range to inner edge of none B 

CZ (RCZi) and range to 

outer edge of CZ (RCZo) none C 

may be displayed by use 
of User Control keys B 
and C respectively. 

Ranges are in nautical 
miles . 



10.0 



RCZi 

RCZi 

RCZo 
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storage Allocation for CZ Range Programs 
Registers ; 



RO; 


°1 


/ EAr 


SO: 


^1 


S: C- / cos 0 

0 rmin 


Rl: 


^1 




SI: 


^2 


B: 6^^^ / 0 

SLD ' rmin 


R2 




^ ^S 


S2 : 


^3 


C : cos 0 

w 


R3: 






S3: 


^4 


D: cos Q 

0W+1 


R4 : 


°3 


/ S 


S4 : 




E: r 

w 


R5: 






S5 : 




I : Control 


R6 ; 




/ r . / RCZi S6: 

^ min ' 






R7 : 






S7 : 






R8 ; 




/ r^ / RCZc 


t S8 : 






R9 : 


^0 




S9 : 






Initial Flag Status 


and Use : 






0 : 


Off, 


Unused 




2: 


Off, set prior to rQ 

calculation in Deep/Deep 
program only. 


1: 


Off, On if RCVR 
in Shal/Shal or 
program only. 


is shallow 
Crosslayer 


3. 


Off, set by data entry 
until r^ is found in 

S/S or Crosslayer program 
only . 



Display Status; DSP 1 



User Control Keys : 



A: 


Data entry 


a : 


B: 


Display 


RCZi 


b • 


C: 


Display 


RCZo 


c : 


D: 






d: 


E: 






e : 
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Step Keys Code 


Explanation 


S tep 


Keys Code 


Explanation 


eg] 


j^LdLh 




21 


11 ; 




051 


2/S 


Si 


Enter DR^^ DS . 


m 


570 1 




35 


45 


Enter data 


052 


k’CLO 


36 00 


. ^ . j 




ISZl 


16 


26 


46’ 




053 


- 


-45 




m 


R'S 






51 




054 


RCL 1 


36 45 




005 


0 






“6£ 


Compute 


055 


V 

RCLl 


-35 


Compute 

O 


006 






00 


Gradients 


056 


36 01 


007 








62 


'So 


057 


+ 


-55 




000 


STo; 




35 


45 


058 


ST02 


35 02 


f 


00? 


ISZl 


i S 


26 


46 




059 


R* 


-31 




010 


RCL3 




36 


63 




060 


RC10 


36 00 




011 


P.CLl 




36 


61 




061 


- 


-45 


Compute C„ 


012 


RCL2 




35 


02 


Si 


062 


RCL ; 


36 45 


iV 


01J 


RCLd 




36 


00 


063 


V 


-35 




014 


(?bOa 


23 


16 


11 




064 


RCLl 


36 01 




015 


RCL5 




36 


05 




065 


4- 


-55 




016 


RCL3 




36 


03 




066 


5704 


35 04 




0ir 


RCL4 




36 


04 


S2 


067 


SF2 


l6 


Initialize 


010 


RCL2 




JO 


02 




068 


1 


61 


AP routine 


019 


GSBn 


23 


lb 


11 




069 


STOC 


35 i3 


020 


RCL7 




JO 


67 




070 $lSl2 


21 02 




021 


RCL5 




36 


05 




071 


GSBc 


23 16 13 




022 


RCLb 




36 


06 


^3 


072 


RCL 3 


36 03 




023 


RCL4 




36 


04 


073 


RCLC 


36 13 


2 A^-i 


024 


CSSq 


23 


16 


11 




074 


eSEol 


23 16 14 


1 


025 


RCLl 




36 


01 




075 


STO0 


35 00 




02b 


RCL7 




36 


07 




076 


ST+0 


35-55 00 




027 


RCLS 




36 


00 


S4 


077 


RCL 3 


36 03 




028 


P.CL6 




36 


06 


078 


EHTt 


-21 




029 


esfia 


23 


16 


11 




079 


EHTf 


-21 




030 


esBk 


23 


16 


12 


Set 1=10 


080 


RCLS 


36 05 


2 A7*2 


031 


Rtid 




36 


00 




081 


RCLO 


36 14 


032 






lb- 


•42 


Mixed layer ? 


082 


GSBcl 


23 16 14 




033 


CT00 




22 


00 




083 


S" + 0 


35-55 00 




034 


ST08 




35 


00 


If no ML, 


084 


ST*0 


35-55 00 




035 


STOl 




22 


01 


0 to R8 


085 


RCLS 


36 05 




03b 


*L8L0 




21 


00 




086 


EHTf 


-21 




037 


RCLS 




36 


09 




087 


EHTf 


-21 




030 


X 






-35 




088 


RCL7 


36 07 


2 A^^ 


039 


CHS 






-22 


= 0 


089 


RCLD 


36 14 


040 


RCLl 




36 


01 


090 


GSB6 


23 16 14 




041 


+ 






-55 




091 


57+0 


35-55 00 




042 


STO0 




35 


11 




092 


57+0 


35-55 00 




043 


QSBc 


23 


16 


13 




093 


RCLl 


36 01 




044 


RCLH 




36 


11 




054 


RCLD 


36 14 




045 


1 






01 


ML present, 


095 


T 


-24 




04b 


CSBcl 


23 


16 


14 


2 Ar-Q to R8 


096 


RCL : 


36 45 


9 A 


047 


STG8 




35 


00 


097 


r 


-24 




048 


ST+S 


35-55 


00 




098 


RCLD 


36 14 




049 


$LBL1 




21 


01 


Set 1=10 


099 


GSBe 


23 16 15 




050 


eSBh 


23 


16 


12 


100 


V 


-35 





CZ Range Program (Deep/Deep Case) 
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step Keys Code 



Explanation 



Step Keys Code 


Explana tion 


151 


OSBoi 


23 16 14 




152 


ST*8 


35-55 08 


to Tg 


153 


DSZI 


16 25 46 


154 


SSBc 


23 16 13 




155 


RC L4 


36 04 1 




156 


RCL8 


36 11 


to r . 


157 


GSBcl 


23 16 14 


min 


158 


ST-6 


35-45 06 




155 


1 


01 




160 


8 


08 


Convert 


161 


5 


05 


RCZi & RCZo 


162 


o 


02 


to nm 


163 


5746 


35-24 06 




164 


5748 


35-24 08 




165 


GSSfc 


25 16 12 


Set 1=10 


166 


♦■L8L6 


21 12 




. 167 


RCL6 


36 06 


Display RCZi 


168 


8/S 


51 




168 


^LBL l 


21 13 


i 


178 


RCLB 


36 08 


Display RCZo 


171 




51 


■ 172 


tLSLa 


21 16 11 




173 


- 


-45 




174 


Ri 


-31 


Gradients 


175 


- 


-45 


Sudroutine 


Id I 1 b 


Rf 


16-31 




ei77 


f 


-24 




178 


570 1* 


35 45 




179 


I5ZI 


16 26 46 




180 


RTS 


24 




181 


fUBLh 


21 lt)“ 12 




182 


1 


di 1 


Set 1=10 


183 


0 


Be 


Subroutine 


184 


STOI 


55 4tj 1 


185 


RTN 


24 




186 


*LBLc 


21 16 13 


Initiation 


187 

188 


PCLl 

EHTt 


36 01 
-21 


189 


EHTt 


-21 


Subroutine 


190 


RTH 


24 




191 


$LBLcl 


21 16 14 




152 


STOC 


35 13 




193 


X 


-35 




194 


X7V 


-41 


Subroutine 


195 


T 


-24 




196 


5700 


35 14 




. 197 


GsBi 


23 16 15 




198 


RCLC 


36 13 




1 99 


C-SSe 


23 16 15 




200 


- 


-45 





lei 

102 



sr+tf 

ST+(3 



35-55 00 
35-55 00 



i03 

"rer- 

i05 

I0€ 



SSBo 



w 

108 

109 

lie 

111 

112 



—FTT 

ST03 

ST04 

tuBLj 

RCL0 

ST*6 

1 

STOE 



23 16 12 

is~2rw 

22 03 
22 04 



Set 1=10 



Pq found ? 



— irsj 

36 00 
35-55 08 
01 
-55 
35 15 



Store Pq 
& first p. 



w 



113 

114 

115 

116 
117 



118 

119 

inr 

121 

122 

123 

124 

125 

126 



*LBL4 

RCL0 

RCLE 

X4V? 

em 



21 94 
36 00 
36 15 
16-35 

22 05 



'min 



found ? 



R4 

STOE 

TTIF 

• 

5 

+ 

ST OB 
COS 
STOC 



-31 
35 15 

-62 
05 
-55 
35 12 
42 
35 13 



cstore 
next Pw 



Increment 
®SLD - 
®3LD 



127 

■nr 

129 

IW 

131 

132 

133 

134 

135 



■ST02 

sro6 



22 02 

irw 

35 06 



Next p 



Store p 



min 



^LB 

5 

COS 

STO0 



-62 
05 
-45 
42 
35 11 



Initiate 

corrections 



136 

137 
136 

139 

140 

141 



OSBc 

RCL2 

1 

GSBcl 

ST+8 

DSZI 



23 16 13 
36 92 
01 

23 16 14 
35-55 08 
16 25 46 



^^3 

corrections 
to Pq 



142 

143 

144 

145 

146 

147 



bS0G 

RCL2 

RCL8 

SSB6 

ST-6 

DSZI 



23 16 13 
36 92 
36 11 
23 16 14 
35-45 06 
16 25 46 



to r 



min 



148 

149 

150 



SSBc 

RCL4 



23 16 13 
36 04 
01 



Ar^ 

corrections 
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SteiD Keys Code 


Explanation 


J tep K eys Code 


Explanation 










201 


X 


-35 








202 


RCLC 


36 13 








203 


7 


-24 


Af 






204 


RCLi 


36 45 


Subroutine 






205 


1S2I 


16 26 46 


(continued) 






206 


T 


-24 








207 




16 31 








208 


RTH 


24 








20S 


$LSLe 


21 16 15 








210 




53 








211 


CHS 


-22 


cos to sin 






212 


I 


01 


Subroutine 






213 


+ 


-55 








214 


7X 


54 








215 


RTH 


24 








216 


R'S 


51 









CZ Range Program (Deep/Peer) Case) 



Step 


Keys Code 


Explanation 


8t*J 




2i 11 




602 

0fiJ 


5 T 0 / 
ISZl 


1 6 


35 45 
2o 46 


Enter data 


m 


R/S 




51 




005 
00^ 
00 r 


5 

C 




-62 

dd 

u2 


Compute 

Gradients 


000 


9T0i 




35 45 


Sn 


009 


ISZI 


16 


26 46 


u 


010 


RCL3 




36 U3 




i 011 


PXil 




36 01 




; 012 


RCL2 




36 02 


Si 


j 0J3 


RCLO 




36 00 


1 


' 014 


QSBa 


23 


16 11 




015 


RCl5 




36 05 




010 


RCL3 




36 03 




dir 


PCL4 




36 04 


oo 


die 


RCL2 




36 92 


C 


015 


QSBa 


23 


16 11 




020 


RCL7 




36 0F 




021 


FCL5 




36 05 




022 


RCLS 




36 06 


Srj 


023 


RCL4 




36 04 


J 


024 


hSEa 




16 11 




025 


RCLl 




36 01 




926 


RCL7 




36 07 




927 


RCLS 




36 06 


§4 


920 


RCL6 




36 96 


929 


5SB<k 


23 


16 11 




939 

931 


Q 

STOI 




09 
35 46 


Set 1=9 


932 


RCLS 




36 00 




933 


RCLS 




36 99 




934 

035 


X 

CHS 




-35 

-22 


Compute Cq 


036 


RCLi 




36 01 




937 


+ 




-55 




938 


STC0 




35 11 




039 


es8c 


23 


16 13 




040 


RCLh 




36 11 




041 

942 


1 

b5B<j/ 


23 


01 
18 14 


2 ^r^ to R8 


043 


STOS 




35 08 




044 


ST+S 


35- 


-55 06 




045 


CSBh 


23 


l8 12 


Set 1=10 


046 


R-' S 




51 


Enter DR DS 


~W7 

048 


x>yo 

xrv 




18-34 

-41 


Reciprocity 

test 


045 

050 


RCL9 




38 00 
-45 


Compute C^ 



Step K eys Code 



Explana tion 



051 


RCLS 


36 05 




052 


V' 


-J5 




053 


RCLl 


36 01 




. 054 




-55 




055 


S702 


35 02 




058 


Ri 


-3l 




957 


RCL9 


36 00 




058 

059 


X>V’-'' 

8700 


16-34 
22 00 


Set Flag 0 


080 


8701 


22 01 


for shallow 


081 


4LBL9 


21 00 


Receiver ■ 


082 


SFl 


16 21 01 




063 


DS2I 


18 25 48 


and 


084 


*LBL1 


21 01 


1 


085 

088 


RuL ! 


-45 
38 45 


Compute 


067 


X 


-35 




96b 


RCLl 


38 01 




085 


* 


-55 




030 


5704 


35 04 




8F1 


QSB6 


23 16 12 


Set 1=10 


972 


*LBL2 


21 02 




973 


GSBc 


Cb 16 lb 




97 4 


RCL3 


36 03 




975 


RCLC 


36 13 


2 Ar. 


976 


6SB6 


23 16 14 


L 


977 


ST 00 


35 09 




970 


S' 7+0 


35-55 00 




. 979 


RCLS 


38 03 




969 


EHT 1 


-21 




. 8S1 


EHTT 


-21 




082 

083 


RCLS 

RCLO 


36 05 
38 14 


2 A772 


084 


GSScI 


23 18 14 




085 


ST+0 


35-55 00 




088 


ST-^9 


35-55 00 




067 


RCLS 


38 05 




90S 


EHTf 


-21 




085 


EHTf 


-21 




050 

051 


RCL7 

RCLC 


36 07 
36 14 


2 


052 


QSB6 


23 18 14 




053 


ST +9 


35-55 00 




054 


ST+0 


35-55 00 




B95 


RCLl 


38 01 




058 


RCLO 


38 14 




853 


+ 


-24 




098 

055 


8CL 1 


36 45 
-24 


2 A 134 


100 


R'ClD 


38 14 
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ir.1 
m * 
A 




Step Keys Code 


Explanation 


181 cs£e 23 16 15 

182 X -35 

183 ST-i-d 35-55 88 

184 5T/-8 35-55 88 




185 0SS6 23 16 12 


Set 1=10 


186 F3? 16 23 83 

187 GT03 22 83 

108 GT04 22 84 


r^ found ? 


109 fLSL3 21 0J 

110 RCL8 36 08 

111 ST+8 35-55 88 

112 1 81 

113 + -55 

114 STOE 35 15 


Store Vq 
2c first r 

w 


115 *LBL4 21 84 

116 RCL8 36 88 

117 RCLE 36 15 

lie 16-35 

119 GT05 22 85 


r^. found ? 
min 


128 R4 -31 

121 STOE 35 15 


Store 
next rv/ 


122 RCLE 36 12 

1 23 . -62 

124 5 85 

125 + -55 

126 ST OB 35 12 

127 COS 42 

128 STOC 35 13 


Increm ent 
®SLD * 
®S1D 


129 GT02 22 82 


Next r rout in 


i38 *1615 21 8S 

131 ST06 35 86 


Store r . 

min 


132 RCLB 36 12 

133 . -62 

134 5 85 

135 - -45 

136 COS 42 

137 STOfi 35 11 


Initiate 

corrections 


138 GSBc 23 16 13 

139 RCL2 36 82 

140 1 81 

141 DSZI l6* 25 46 

142 OSBal 23 16 14 

143 DSZI 16 25 46 


; 

corrections i 

j 

to Tq 1 


144 ST +8 35-55 88 


145 gSBc 23 16 13 

146 RCL2 36 82 

147 RCLh 36 11 

148 GSBd 23 16 14 

149 Sr+b 55-55 86 


to r . 
min 


158 FI? 16 23 81 





Step K eye: Code 



151 

J52 

J55 

154 

155 

156 

157 
15$ 
15.9 



05ZJ 

esec 

RCL4 

1 

GSBcl 

OSZI 

FI? 

CHS 

ST+8 



16 25 46 
23 16 13 
36 $4 
01 

23 It) 14 
l6 Z5 46 
16 23 81 
~22 
35-55 88 



168 

161 

162 

163 

164 

165 

166 



GSBc 

RCL4 

RCL8 

Gi'Bd 

FI? 

CHS 

ST-6 



23 16 13 
36 84 
36 1 1 
23 16 14 
16 23 81 
-22 
55-45 86 



Jbr 

lb& 

169 

178 

171 

172 



1 

S 

5 

2 

ST ^6 



81 

88 

85 

82 

35-24 86 
35-24 06 



173 tS61) 23 l6 12 



1/4 

175 

176 



*lblb 

RCL6 

R/S 



21 12 
36 86 
51 



TT 
1 7$ 
el 79 



$LBLC 

RCLS 

pc 



21 13 
36 88 
5i_ 



Ex plana ti on 



corrections 
to r„ 



to r 



min 



Convert 
RCZi £c RCZo 
to nm 



Set 1=10 



Display RCZi 



Display RCZo 



186 *LBu 21 i6 11 

1$1 - -45 

1S2 Rl -31 

183 - -45 

184 Rf 16-31 

185 ^ -24 

1S6 STOi 35 45 

187 I5Z1 16 26 46 

183 RTN 24 



Gradients 
Subrout ine 



189 

198 

191 

192 

193 



*LBLi> 

1 

8 

STOI 

RTN 



21 16 12 
81 
88 
J5 4t> 
24 



JP4 


*LBLc 


21 16 13 


195 


RCLl 


36 81 


IPb 


EfiTf 


-21 


197 


EHTf 


-21 


198 


RTN 


24 


J QQ 


rLBLd 


21 16 14 


280 


sroc 


35 13 



Set 1=10 
Subroutine 



Initiation 

Subroutine 



CZ Range Program (Shal/Shal <3: Crosslayer Cases') 



90 




91 



User Instructions for CZ Gain and Transmission Loss Program 



Step 



Instructions 



Input Keys 



1 



2 

3 



4 

5 



6 



Run one of the CZ Range 
Programs . Leave the 
calculator on and all of 
the data storage regis- 
ters unchaged. 

Load the CZ Gain and 
Transmission Loss Program. 



Enter a zero if the Shal/ 0 A 

Shal or Crosslayer program 

was used, or enter a one 

if the Deep/Deep program 1 

was used. 



Enter a frequency of in- f (Hz) B 

terest in Hz. 



Enter the number of a CZ n C 

annulus of interest. (Step 
5 may be repeated for as 
many CZ annuli as desired.) 



To determine the TL 

n 

values for different 
frequencies , return to 
step 4 . 



Output 



G (dB) 



a (dB/m) 

TL (dB) 
n 
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storage Allocation for CZ Gain and Transmission Loss Program 



Registers : 



RO: 


Ar 


SO: 


^1 


A: 


®S0 ^ ®R0 


Rl: 


^1 


SI: 


^2 


B; 


®SLD ^ ®rswp 


R2; 




S2; 


53 


C; 


cos 0 / 

w ' 


R3; 


^2 


S3 : 


54 




cos 9 

rswp 


R4 ; 




S4 : 


a terms 


D: 


cos 9 , T 
w+1 


R5 : 




S5: 




E: 


^0 ^ ®SR ^ ® 


R6 ; 


RCZi 


S6 : 




I; 


Control 


R7 ; 


C4 / a 


SI: 








R8 ; 


RCZO 


S8: 








R9 : 


G terms / G 


S9 ; 









Initial Flag Status and Use ; 



0; Off, On for shallow 
sound source. 


2: 


Off, unused. 


1: Off, On to decrease the 

increments used in 
finding . 


3: 


Off, set by data entry 
until r^ found. 


Display Status; DSP 0 






User Control Keys: 






A; Compute G 


a : 




B : Compute a 


b : 




C; Compute TL^ 


c : 





Steo Keys Code 


explanation 


S tet> 1 ' eyr. 1 o < 


Le 


iXX[.lana t. iof 










001 

002 

663 


*(.Slh 

.V=0'- 

SF6 


Jj il 
it? 21 6tf 


3 


et FO for 
hal Source 


051 

052 
05J 


i 

+ 

Q T 03 


22 


01 

-55 

63 


Large © 3 ^ 
increment 


004 


1 


61 


Initialize 


054 


tiSL 1 


21 


01 




005 


STOC 


35 13 


r 


routine 


055 


SFl 


16 21 


01 




000 


*lBL6 


21 <30 






05o 


RCLB 


36 


12 


first small 


007 


FCLl 


Jo 0] 






057 






*62 


increm ent 


000 


EtiJi 


-21 






050 


5 




09 


^SLD 


009 


EHTt 


-21 


2 


Ari 


059 


- 




-45 




616 


RCL3 


36 63 






0o‘0 


OT03 




03 




611 


RCLC 


36 13 






661 


*LBL2 


21 


62 




612 


QSBa 


23 16 11 






062 


;v>'r? 


lo- 


■34 


®rswp • 


613 


ST 06 


35 06 






663 


erc4 




64 




614 


RCL3 


36 63 






664 


RCLB 


36 


12 


3ma.'Ll ft 


615 


EHTt 


-21 






665 


« 




■62 


w ^ T* 


616 


EHTt 


-21 






066 


1 




61 


increment 


617 


RCL5 


36 65 


2 




667 






-55 




616 


RCLD 


36 14 




668 


«LBL3 


■ 21 


0J 




619 


6S8a 


23 16 11 






669 


STOB 


J5 


12 


Store nev; 


626 


ST+6 


35-55 06 






676 


COS 




42 


®SLD ^ 


621 


RCL5 


36 65 






671 


STOC 


35 


1 J 


COS G .,-, 


622 


EHTt 


-21 






672 


BTO0 


22 


00 


oiD 


623 


ENT1 


-21 






673 


*L3l4 


21 


04 




624 


RCL7 


36 67 


2 




674 


1 




01 




625 


RCLD 


36 14 






675 


0 




00 


Convert 


626 


eSB<x 


23 16 11 






676 


5 




05 


RCZi & RCZo 


627 


ST-i-6 


35-55 66 






677 






02 


to meters 


628 


RCLl 


36 01 






678 


srx 6 ' 


35-35 


0o' 




629 


RCLD 


Jo 14 






679 


STxS 


35-35 


00 




636 


T 


-24 






686 


RCLB 


36 


12 


^ "6 0 T' 0 


631 


RCL! 


36 45 






681 


COS 




42 


0 U U i. c 


632 


r 


-24 


2 


Ari. 


682 


STOC 


35 


IJ 


cos G 

rswp 


633 


RCLD 


36 14 






683 


PCL2 


36 


02 




634 


QSBh 


23 16 12 






684 


X 




-J5 




635 


X 


-35 






685 


RCLl 


36 


01 


1 G 


' 636 


ST+0 


35-55 66 






686 


T 




■24 


^3R 

1 


637 


ST+0 


35-55 86 






087 


COS- 


16 


42 


1 

( 

1 


638 


1 


61 




6SS 


STOE 


35 


15 


I 



S39 

040 



641 

642 

643 



644 

645 

646 

647 



646 

645 

W 



a 

ST 01 



00 
J5 46 



RCL6 

F37 

STOE 



36 00 
16 23 63 
35 15 



RCLE 

XZY 

FI'-' 

0TC2 



36 15 
-41 
16 23 61 
22 62 



t / UJ 



16-34 
22 61 



PC LB 



JO 12 



•Set 1=10 



Store r 



0 



Checic for 
Small inc . 

®SLD 



Decrease 
©SLD inc . 



889 


RCLE 


36 02 


690 


RCLl 


Jo 01 


091 


r 


-24 


692 


COS" 


16 42 


693 


Fd'^ 


16 23 00 


694 


CHS 


-22 


095 


STOP 


J5 11 


696 




-55 


697 


0 -K’ 


lo' 45 


69S 


ST 09 


J5 09 



oO 

(negative if 
source shal) 



655 

100 



RCLE 

PlLh 



36 15 
36 11 



A© 



CZ Gain and Transmission Loss Program 



94 



10] 


- 


-45 




10E 




02 




103 


T 


-24 


COS 


104 


COS 


42 




105 


srx5 


35-35 09 




106 


RCLC 


36 13 




107 


RCL4 


36 04 




108 


X 


-35 


Q 


109 


RCLl 


36 01 


^RR 


110 


T 


-24 




111 


COS- 


16 42 




112 


RCL4 


36 04 




113 


RCLl 


36 01 


Q 


114 




-24 


®RO 


115 


cos-‘ 


16 42 





step Keys Code 



Explanation 



Step Keys Code 


Explana tion 


151 


4 


04 




152 


1 


01 




153 


e 


00 


Second term 


1 54 


ft 


60 


155 


PCIA 


36 Ji 




156 


-f 


-55 




157 


•r 


-24 




158 




56 




159 


RCLl 


36 56 


Store a 


160 


STG7 


35 07 


161 




16 56 


(Clear S4) 


162 


1 


01 




163 








164 


d 


66 


Convert a 


165 


5 


6S 


from dB/kyd 


1 66 


4 


04 


to dB/m 


1 6 f' 


EEX 


-23 




168 


J 


63 




169 


CHS 


—*141 




170 


STxr 


35-35 6f 




171 

172 


RCL7 
R - S 


36 97 
51 


Display a 


173 


*LBLC 


21 13 




174 

175 


RCLS 

X 


36 06 
-35 


nRCZi 


176 


EHJt 


-2l 




177 


log 


16 32 




178 




02 


201og(nRCZi) 


179 


0 


00 


130 




-35 




Ibl 


•-/ 

i 


-41 




1S2 


RCL7 


3o 97 


_a (nRCZi) 


183 


X 


-35 ' 




18^ 




-55 




1S5 


RCLS 


36 09 


Compute & 


186 


- 


-45 


Display TL^ 


187 


Ry’S 


51 


16S 


tLBLci 


21 16 11 




169 


STCC 


35 13 




190 


\ 


-35 




191 


X^r 


-41 




192 




-24 




193 


sroo 


35 14 


Ar 


194 


QSSii? 


23 16 12 


Subroutine 


195 


RCLC 


36 13 




196 


QSBh 


23 16 12 




197 


- 


-45 




198 


V 


-35 




1 55 


PC Li 


36 13 




20P 


T 


-24 





lib' 

Jir 

ns 



-55 

02 

-24 



112 

120 



SIH 
ST ^9 



41 

35-24 09 



sin 0, 



121 RCLS 

122 RCLS 

123 STx9 



36 0S 
36 06 
35-35 09 



RCZi 



124 

125 



ST^9 



-45 
35-24 09 



czv; 



126 

127 

128 

129 

130 

131 

132 



RCL9 

LOG 

1 

0 

X 

ST 09 
RyS 



36 09 
16 32 
01 
00 
-35 
35 09 
51 



Compute and 
Display G 



133 *LBLB 

134 EEX 

135 3 

136 4 



21 



12 

-23 

03 

-24 



Convert f 
to IcKz 



137 

138 



STOrt 



53 
35 11 



Store 



X 



139 




-62 


140 


i 


01 


141 


X 


-35 


142 


I 


01 


143 


/?clh 


36 11 


144 


+ 


-55 


145 


T 


-24 


146 




56 


147 


4 


64 


148 


0 


66 


149 


ECLh 


36 1 i 


150 


X 


-35 



Attenuation 

Coefficient 



First term 
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Step Keys Code 


Explana tion 


Step Keys ^o(]e 


Sxp is ri;j 1 ion 


201 RCLi 36 45 

202 ISZI 16 26 46 

203 - -24 

204 hBS 16 31 

205 2 02 

206 X -35 

207 RTti 24 


Subroutine 
(Continued) 1 






206 *LBL6 21 16 U 

209 X'2 55 

210 CHS -22 

211 1 01 

212 + -55 

215 -TX 54 

214 RTN 24 

215 RyS 51 


cos to sin 
Subroutine 
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